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ABSTRACT 

Recent numerical simulations have suggested the possibility of forming double intermediate mass 
black holes (IMBHs) via the collisional runaway scenario in young dense star clusters. The two IMBHs 
formed would exchange into a common binary shortly after their birth, and quickly inspiral and merge. 
Since space-borne gravitational wave (GW) observatories such as LISA will be able to see the late 
phases of their inspiral out to several Gpc, and LIGO will be able to see the merger and ringdown out 
to similar distances, they represent potentially significant GW sources. In this Letter we estimate the 
rate at which LISA and LIGO will see their inspiral and merger in young star clusters, and discuss 
the information that can be extracted from the observations. We find that LISA will likely see tens 
of IMBH-IMBH inspirals per year, while advanced LIGO could see ~ 10 merger and ringdown events 
per year, with both rates strongly dependent on the distribution of cluster masses and densities. 
Subject headings: stellar dynamics — black hole physics — gravitational waves 



1. INTRODUCTION 

Observations suggesting the existence of intermediate- 
mass black holes (IMBHs) have mounted in recent 
years. Ultra-luminous X-ray sources (ULXs) — point X- 
ray sources with inferred luminosities > 10 39 erg/s — may 
be explained by sub-Eddington accretion onto BHs more 
massive than the maximum of ~ lOM© expected from 
stellar core collapse (Miller & Colbert 2004). Similarly, 
the cuspy velocity dispersion profiles in the centers of 
the globular clusters M15 and Gl may also be explained 
by the dynamical influence of a central IMBH (van der 
Marel et al. 2002; Gerssen et al. 2002; Gebhardt et al. 
2005), although this conclusion remains somewhat con- 
troversial (Baumgardt et al. 2003). 

The most likely formation scenario for an IMBH is the 
collapse of a very massive star (VMS) , which was formed 
early in the lifetime of a young star cluster via a runaway 
sequence of physical collisions of massive main-sequence 
stars (Portegies Zwart et al. 1999; Ebisuzaki et al. 2001; 
Portegies Zwart & McMillan 2002; Giirkan et al. 2004). 
This scenario has been studied in detail for star clus- 
ters without primordial binaries, with recent work show- 
ing that runaway growth of a VMS to ~ 10 3 Mq oc- 
curs generically in clusters with deep core collapse times 
shorter than the ~ 3 Myr main-sequence lifetime of the 
most massive stars (Freitag et al. 2006). 

Due to the computational cost of simulating the more 
realistic case of star clusters with primordial binaries, 
it is only recently that such simulations have been per- 
formed (Portegies Zwart et al. 2004; Giirkan et al. 2006). 
The work of Giirkan et al. (2006) was the first to system- 
atically study the influence of primordial binaries on the 
runaway growth process. They showed that stellar colli- 
sions during binary scattering interactions offer an alter- 
nate channel for runaway growth, with the main result 
that clusters with binary fractions larger than » 10% 
generically produce two VMSs, provided the cluster is 



sufficiently dense and/or centrally concentrated to trig- 
ger the runaway earlier than ~ 3 Myr in the absence of 
primordial binaries. Observations and recent numerical 
calculations suggest that star clusters may be born with 
large binary fractions (> 30%; Hut et al. 1992; Ivanova 
et al. 2005), implying that all sufficiently dense and mas- 
sive star clusters could form multiple VMSs. 

The VMSs formed will undergo core-collapse su- 
pernovae and likely become IMBHs on a timescale 
of ~ 4 Myr after cluster formation (the lifetime of a 
VMS is extended slightly by collisional rejuvenation; 
see, e.g. Freitag et al. 2006). After their separate 
formation, the two IMBHs will quickly exchange into 
a common binary via dynamical interactions. The 
IMBH-IMBH binary (IMBHB) will then shrink via 
dynamical friction due to the cluster stars, on a 
timescale ~ t r (m) /Mimbh IS 10 Myr, where t r is the 
core relaxation time, (m) is the local average stellar 
mass, and (m) /A/jmbh < 10~ 2 . Note that since 
t r scales inversely with (m) for fixed core velocity 
dispersion and mass density, the dynamical friction 
timescale is independent of (m) (Binney & Tremaine 
1987). The IMBHB will then shrink further via dy- 
namical encounters with cluster stars (Quinlan 1996; 
Yu & Tremaine 2003; Miller 2005), until it merges 
quickly via gravitational radiation, on a timescale 
1 Myr (cr c /20kms- 1 ) 3 (p c /10 5 M pc^-^Mtmbh/IO 3 M q )-\ 
where a c is the cluster core velocity dispersion and p c 
is the core mass density (Quinlan 1996, eqs. [29] and 
[30]). This timescale has also been confirmed by 
numerical scattering calculations (Gultekin, private 
communication) . 

Only the more massive IMBHBs merge in the LISA 
band of 10 _4 -1 Hz (redshifted binary mass M z = (1 + 
z)M >4x 1O 3 M , where M is the total binary mass). 
Fig. 1 shows the final gravitational wave (GW) frequency 
// (the frequency at the inner-most stable circular orbit 
if within the LISA frequency range (large M z ), otherwise 
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Fig. 1. — The final GW frequency ff (see text), and the frequency 
1 yr prior, fi, for an IMBHB with total mass M and reduced 
mass parameter 77, as a function of redshifted binary mass M z , for 
r) = 0.25 (equal-mass binary) and 77 = 0.1 (mass ratio 0.13). (The 
final frequency is roughly independent of r\.) 
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Fig. 2. — Luminosity distance, to which an IMBHB of 

total mass M and reduced mass parameter 77 can be seen via its 
inspiral with LISA with S/N = 10 for a 1 yr integration, and via 
its merger and ringdown with S/N = 8 for iLIGO and AdLIGO, as 
a function of the redshifted mass M z . The corresponding redshift 
(calculated using the WMAP year 3 cosmological parameters, as 
discussed in the text) is shown on the right vertical axis. 



the maximum LISA frequency of w 1 Hz (small M z ), as 
in Will (2004)), and the frequency 1 yr prior, fi, as a 
function of redshifted mass M z , for the reduced mass pa- 
rameters r] = 0.25 (equal-mass binary) and r\ = 0.1 (mass 
ratio 0.13) (see, e.g., Will 2004). For a wide range in M z , 
the late stages of inspiral clearly span the LISA "sweet 
spot" (roughly a decade centered on 10~ 2 ' 2 Hz), implying 
that LISA could easily detect the chirp signal, enabling a 
measurement of the masses of the binary members. Such 
an observation would be direct evidence for an IMBH. 

In Sec. 2 we estimate the rate at which LISA will ob- 
serve inspiral of IMBHBs in young star clusters. In Sec. 3 
we estimate the rate at which LIGO will observe their 
merger and ringdown. Finally, in Sec. 4 we discuss the 
observational consequences. 

2. ESTIMATING THE LISA DETECTION RATE 



We first need to know the distance to which LISA can 
see IMBHB inspirals. Following the techniques in Will 
(2004) and Flanagan & Hughes (1998), we adopt the lat- 
est LISA sensitivity curve (Larson 2003), including con- 
fusion noise from Galactic white dwarf binaries (Bender 
& Hils 1997), and calculate the maximum luminosity dis- 
tance, <1l(z), to which an IMBHB of total mass M and 
reduced mass parameter r/ can be seen with S/N =10 
for a 1 yr integration. The results are shown in Fig. 2 as 
a function of M z , for r\ = 0.25 and r/ = 0.1. Note that 
the results of Giirkan et al. (2006) show that the masses 
of the IMBHs never differ by more than a factor of a 
few (rj > 0.15). Thus LISA will be able to see typical 
IMBHBs (M - 10 3 M Q ) out to a few Gpc. 

With this information in hand, we first make a crude 
estimate of the LISA event rate. Following Miller (2002), 
we write for the total rate 
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The first factor, J Zmax (dV c / dz)dz , is the integrated co- 
moving volume of space in which LISA is sensitive to the 
events. The second factor, dN c \/dV, is the number den- 
sity of star clusters sufficiently massive to form IMBHBs. 
Since the globular clusters we currently see were likely at 
least a few times more massive at formation (Joshi et al. 
2001), we set this factor to the current density of globu- 
lar clusters in the local universe, dN c \/dV w 8/i 3 Mpc -3 
(Portegies Zwart & McMillan 2000). The third factor, g, 
is the fraction of sufficiently massive clusters that have 
a large enough binary fraction and initial central den- 
sity to produce IMBHBs. Since initial cluster structural 
parameters are largely unknown, we treat g as a param- 
eter. The fourth factor is the event rate per IMBHB- 
producing cluster, taken to be one divided by the age of 
the universe, since only one IMBHB is formed per clus- 
ter over its lifetime. We adopt a ACDM cosmology, with 
parameters Qm = 0.24, f^A = 0.76, and h = 0.73, for 
which tu = 13.8 Gyr (Spergel et al. 2006). Putting this 
together for d^ = 4.9 Gpc (z max = 0.79), the distance 
to which LISA can see IMBHBs with M = 2 x 10 3 M Q , 
cq. (1) gives Rtv 1(5/0.1) yr" 1 . 

Writing down a generalized form of the rate integral in 
eq. (1) is straight-forward. Since the time between cluster 
formation and IMBHB merger is we assume that 

the merger is coincident with cluster formation. Thus 
the rate integral is 



R = 
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Here R = dN cvcnt / dt is the event rate observed at z = 
by LISA, d 2 MsF/dV c dt e is the star formation rate (SFR) 
in mass per unit of comoving volume per unit of local 
time, g c \ is the fraction of star forming mass that goes 
into star clusters more massive than 10 3 ' 5 M© (generally 
a function of z), g is as above, and d 2 N c \/dMsY,c\dM c \ 
is the distribution function of clusters over individual 
cluster mass M c \ and total star forming mass in clus- 
ters Msf.cI- Finally, dt e /dt is simply (1 + z) _1 , and 
dV c /dz is the rate of change of comoving volume with 
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rcdshift, which is a function of cosmological parameters 
(Hogg 1999). Note that we set z max = 5, since this is 
roughly the limit to which the cosmic SFR can be traced. 
Thus the integral in eq. (2) should be considered a mild 
lower limit to the true rate. We now discuss each element 
in eq. (2) in more detail. 

Following Porciani & Madau (2001), we adopt three 
different choices for the SFR: 

/ d 2 M \ 

{d\rdt) SF . = C i h ii6 F{z)G i {z)M Q yr- 1 Mpc" 3 , (3) 

with i = 1,2,3 denoting the different rates, Ci a constant, 
Gi(z) a function of z, he^ = h/0.65, and F(z) = [T2m(1 + 
zf + n k (l + zf + fi A ] 1 / 2 /(l + zf' 2 . The first is from 
Madau & Pozzetti (2000), with d = 0.3 and Gi{z) = 
e 3 - 4z /(e 3 8z + 45), which peaks between z = 1 and 2 and 
decreases at larger redshift. The second is from Stcidel 
et al. (1999), with C 2 = 0.15 and G 2 (z) = e 3Az /(e 3Az + 
22), which is roughly constant for z > 2. The third is 
from Blain et al. (1999), with C 3 = 0.2 and G 3 (z) = 
e 3.05z-o.4^ e 2.93z _|_ -^g^ which increases above z ~ 2. 

Measuring the fraction of star-forming mass in clusters 
is difficult for anywhere but the local universe. Similarly, 
while we know reasonably well the initial cluster condi- 
tions required to form an IMBHB (Giirkan et al. 2006), 
we know much less well the distribution of cluster proper- 
ties at birth. We therefore treat g c \ and g as parameters, 
taking g c \ = 0.1 and g = 0.1 somewhat arbitrarily as 
canonical values. 

Assuming that the spectrum of cluster masses is nei- 
ther a function of cosmic epoch nor the total star forming 
mass available for clusters, the factor d 2 N c \/dM§v^\dM c \ 
can be separated as 

d 2 N cl /(M cl ) 
dM SFM dM cl f M cl f(M cl )dM cl ' {) 

where f(M c \) is the (normalized) distribution function 
of cluster masses. For this we adopt the power-law form 
observed for young star clusters in the Antennae, which 
is thought to be univers al: /(M cl ) oc M^ 2 (Zhang & Fall 
1999), with a lower limit of 10 3 ' 5 Mq and an upper limit 
ofl0 7 M Q . 

It is the limit M c i im ; n (z) in eq. (2) that encodes all 
information about the detectability of an IMBHB inspi- 
ral by LISA. Specifically, the redshift to which LISA can 
see the inspiral is a function of the binary mass, which 
is itself a function of the host cluster mass. Adopting 
an efficiency factor /gc for the fraction of cluster mass 
going into the IMBHB, this relationship is inverted to 
obtain Mci. m i n (z). Recent numerical work shows that 
the efficiency factor is /gc ~ 2 x 10~ 3 , independent of 
cluster initial conditions (Giirkan et al. 2004), which we 
take as our canonical value. At low redshift, M C L m i n (z) is 
clamped at the value M c \ = 200M©//gc, set by adopt- 
ing the definition that an IMBH have mass > 10 2 Mq. At 
high redshift (z > 5, so not relevant to our calculation), 
M c \ tin i n (z) is clamped at the value of 10 Mq from the 
cluster mass function; in other words no cluster is suf- 
ficiently massive to produce an IMBHB massive enough 
to be observable by LISA, so the integral is zero. 

We numerically integrated eq. (2) for the different 
SFRs in eqs. (3), for S/N = 10 and an integration time 
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Fig. 3. — Integrand of the rate integral in eq. (2) for the three 
different SFRs in eqs. (3), for r\ = 0.25 and r\ = 0.1. 



of 1 yr, to find that the rate is 

R(r,= 0.25) « 40-50 (M) ( X) yr-i , ( 5 ) 

with the spread in the coefficient from the different SFRs. 
The coefficient decreases to 20-25 for rj = 0.1. The rate 
is dominated by clusters in the mass range 10 6 -10 6 5 Mq 
(IMBHB mass 2 x 10 3 -6 x 10 3 Mq), with more than half 
the contribution to the rate coming from this mass range, 
for both 77 = 0.1 and 77 = 0.25, and for all three SFRs in 
eq. (3) (except SF3 for 77 = 0.1). Note that eq. (2) is only 
strictly valid when the source is visible by the instrument 
for less than the integration time. This turns out not 
to be precisely correct. A typical IMBHB with mass 
M = JgcIO 6 ' 25 Mq takes roughly 4 years to cross the 
LISA band from the edge of the white dwarf confusion 
knee at ~ 2 mHz to the upper edge of the band at w 1 Hz. 
Thus the rate presented in eq. (5) is an underestimate by 
of order a factor of a few. 

Fig. 3 shows the integrand of the rate integral in eq. (2) 
for the three different SFRs in eqs. (3), for 77 = 0.25 and 
77 = 0.1. Most events originate from z <~ 1. Unfor- 
tunately, neither R nor dR/dz is particularly sensitive 
to the cosmic SFR, with dR/dz decreasing quickly above 
z 2 even when the SFR is increasing (as in SF3). Thus 
observations of IMBHB inspirals will not be very infor- 
mative about the cosmic SFR. However, they will likely 
yield a handle on the fraction of star formation that is in 
compact massive clusters. 

3. ESTIMATING THE LIGO DETECTION RATE 

Shortly after the two IMBHs merge, the merger prod- 
uct can be well described as a single perturbed black hole, 
emitting GWs at its quasinormal frequencies. Largely 
falling within the initial and advanced LIGO (iLIGO 
and AdLIGO) sensitivity bands, the merger and ring- 
down waves will likely carry a few percent of the rest 
mass energy of the hole (see, e.g. Flanagan & Hughes 
1998). Numerical simulations suggest that a merging 
pair of nonspinning equal-mass black holes will emit a 
fraction e ~ 0.03 of their rest mass in merger and ring- 
down GWs, forming a black hole with spin parameter 
a ~ 0.7 (Baker et al. 2002; Campanelli et al. 2006; Baker 
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et al. 2006). Under these conditions, the ringdown fre- 
quency is given by (see eq. [3.17] of Flanagan & Hughes 
1998) 

^2^< 1 - « 1 -°) 3/1 ^ 18 »(T^)" H - 

(6) 

We can express the distance to which we are sensitive to 
ringdown waves at signal-to-noise ratio p as 

where S(f) is the spectral noise density of LIGO. Com- 
bining this expression with the concordance cosmolog- 
ical model and iLIGO and AdLIGO sensitivity curves, 
we find the range to which LIGO can detect ringdown 
shown in Fig. 2. 

To obtain a conservative estimate for the rate at which 
iLIGO and AdLIGO could detect these mergers with 
a ringdown-only search, we use eq. (1) with a moder- 
ately optimistic range of cIl ~ 100 Mpc for iLIGO and 
di = 2Gpc for AdLIGO. The expected detection rate 
is then 10- 4 (.g/0.1)yr- 1 and 1( 5 /0.1) yr~\ for iLIGO 
and AdLIGO, respectively. More detailed estimates us- 
ing machinery analogous to eq. (2) increase these esti- 
mates by roughly an order of magnitude, making the 
rate for AdLIGO 10( 5c i/0.1)(g/0.1) yr" 1 . 

4. DISCUSSION 

It appears likely that LISA will see tens of IMBHB 
inspiral events per year, while AdLIGO could see ~ 10 
merger and ringdown events per year, with both rates 
strongly dependent on the distribution of cluster masses 
and densities. Detection of an IMBHB would have pro- 
found implications. A match-filtered observation of the 



inspiral would yield the redshifted masses of the black 
holes, directly confirming the existence of IMBHs. It 
would also yield the luminosity distance to the source; 
with enough observations, constraints could be placed 
on the cosmic history of star formation in dense, mas- 
sive clusters. Detection of the ringdown signal from the 
merger product will additionally yield its spin, which 
may provide insight into its formation history. 

Typical IMBHBs spend > 10 6 yr inspiraling through 
the LISA band, with nearly all of that time spent at low 
frequencies (< 10~' ! Hz). In the low frequency region 
they will thus appear as a large number of monochro- 
matic sources, possibly contributing to confusion noise 
and increasing the noise floor (e.g., Farmer & Phinncy 
2003). A detailed calculation of this is beyond the scope 
of this Letter. However, we note that if their contribu- 
tion is similar in magnitude to that of Galactic compact 
object binaries (Bender & Hils 1997), the rates predicted 
in eq. (5) would decrease by about 20%. 
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